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, We calculate the CP-violating effects in the top quark semi-leptonic three body decays 

CN ' induced by the super symmetric CP-odd phase of the top squark trilinear soft breaking term 

arg(Aj). The light top squark mass is assumed to be close to the top quark mass mi ~ nit. 
. The CP-conserving phase is provided by the x + an d X° cu t- We find that the partial rate 

asymmetry is in the 0.1% level. In the most favorable parameter region the decay rate 
^ ■ asymmetry can reach up to 0.55%. 

a: 
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^ 1 Introduction 

> : 

ix! ■ 

Top quark physics is sensitive to new physics, which may exist near the electro-weak scale, due 
to its large mass. Experimental and theoretical research about CP violation in top sector is one 
way to reveal new physics. To study the top quark CP-odd effects has its own advantage that 
the uncertainties coming from hardron matrix elements can be avoided. 

If m t is near to m newp hy S ic S , the CP-assymmetry effects in top quark decays can be induced by 
new particles. Until now a lot of works about CP-asymmetry effects in the top quark decays jl] have 
been done within the supersymmetric model. Most of those works have made the assumption that 
the mass of the light top squark is much smaller than that of the top quark. B.Grzadkowski and 
W.Y.Keung calculated the CP violating effects induced by t6g-loop||. This contribution requires 
the condition mi + nig < m t , which has already been excluded. E.Christova and M.Fabbrichesi 
computed the effects induced by £6x°-1° pII which requires m^ + m x o < m t . S. Bar-Shalom et al 
gave the CP asymmetry in top quark decay induced by tx + X°- It is at best 0.3% when light stop 
mass is between ~ 50GeV and ~ 70CreVfl. 
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However, if the light top squark mass mi x is approximately as heavy as the top quark, CP- 
asymmetry effect in top quark two body decay induced by supersymmetric CP-odd phase will not 
be observable because the top squark can not run on shell to produce the necessary absorptive cut. 
In this work we considered this case under the assumption that the light chargino is much lighter 
than rrif . Under this condition, the x°{ which is always assumed to be the LSP ) and x + can 
provide necessary absorptive cut in top quark three body decays, such as in the process t — > bv T r 
considered in the work. These two particles can be on shell in the top quark three body decay 
loop diagrams when the invariant mass of the lepton pair is sufficiently large. To our knowledge, 
a study on CP asymmetry in the top quark three body decays is missed in literature. 

In the present work the mass of the light top squark is assumed to be above 140GeV. Taking 
into account the direct experimental limit on super particles || and the indirect limit coming from 
neutron EDM limit 0, we take \x to be real and , the lightest neutralino to be above 30GeV and the 
lighter chargino to be above 65GeV. The large mass of stop leads to relatively small CP violating 
effects. Nevertheless, the CP-odd effects can reach up to 0.55% in the most favorable parameter 
space. 

The paper is organized as following: in Sec. 2 we analyze the possible new CP violating sources 
in MSSM and present our simplifying assumptions in calculation. In sec. 3 we sketch the main 
steps of our calculations. In sec. 4 we present our numerical results and in sec. 5 and sec. 6 
we discuss and summarize our results. The mass matrices for charginos, neutralinos and squarks 
are given in appendix A. In appendix B we give the relevant pieces of the Lagrangian for our 
calculations and some analytic results are presented in appendix C. 



2 CP violating phases in the low energy supersymmtric 
model 

The most general form of the low energy Lagrangian of MSSM|7|, ||, which is SU(3) x SU (2) x U(l) 
gauge invariant and does not violate the SM conservation laws, can be written as 



where e\2 =-1. The hat ' " ' indicates that the corresponding letter represents a superfield. The 
capital indices I, J denote generations. i,j refer to the components of a SU(2) doublet. The l,d,u 
are the Yukawa coupling matrices. The soft breaking terms can be divided into three pieces, 




(1) 



The superpotential W is given by 



W = eyinHlH? + l IJ E\l]R J + u IJ H*Q I j U J + d IJ HlQ I j D J ) 



(2) 




gaugmo 



+ c 



'trilinear 



(3) 
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These are the scalar particle mass terms, gaugino mass terms and the trilinear soft breaking terms 
respectively. They are given by 



Scalar = 6tfA*.ltffl* - m% Jltf | 2 - mfrjfl? | 2 ~ {m\) IJ L? • Lf - 

{ml) 1 ^ 1 * ■ R J - (mj)"gf ■ Q- - 

(mlY'D 1 * -D J -{mlY'U 1 * ■ U J , (4) 
^g aug i n o = ^(miA B A B + m 2 X i w X i w + m 3 X G X G ) + H.C. , (5) 

Arilinear = £y( (AeI) 1 J H} Lj R J + {A D d) IJ U\Q^D J + 

{A u u) IJ H 2 Q 1 j U I )+H.C. , (6) 

where 7711,7712,771,3 are U(l), SU(2) and SU(3) gauginos masses respectively. Fields in £ sca iar and 
^trilinear are scalar components of the corresponding superfields. 

In general, all the coupling parameters in the above expressions except those of the diagonal 
terms in £ S caiar may be complex which may be the CP violating sources. However, not all of 
them are physical and, even the physical parameters are too many to be disposed of. In actual 
calculations, simplifying assumptions must be made. We get the physical CP violating phases by 
the following steps. 

First, we take the GUT assumption that the m^s are universal at the GUT scale and can 
be set real by a phase rotation||. Thus, the m^s are real at any scale. Second, we adjust the 
global phase between the two Higgs superfields so that /x s is real. This adjustment makes the 
two vacuum expectation values of the neutral Higgs fields v±, v 2 real||. After the adjustment 
the phases of the two Higgs superfields are fixed and \x is complex in general. Third, l IJ , d IJ , 
and u IJ in the superpotential are diagonalized and the unphysical phases are absorbed by quark 
superfields similar to that done in the standard model. This leaves a CP violating phases 5km m 
the Kinetic terms after the superfield are redefined. Fourth, to suppress the FCNC process in the 
SUSY extension of SM, as an approximation we require that all the matrices in £ so ft, m 2 L , m%, 
rrin, mfj, m 2 D , Ajj, Ad and Ae are flavor diagonal in the basis where l IJ , d IJ and u IJ are diagonal 



(flavor alignment ||10||). Then the Hermitian matrices in £ sca i ar are now all real . The phases of all 
squarks are fixed after the third step. Therefore, Ap, Ajj are generally complex. In conclusion, 
5km, arg(/x), a.rg(A D )s and arg(Ajy)s are the CP violating phases in the low energy supersymmtric 
model under our assumptions. 

arg(/i) can not be larger than the order ~ O(10~ 2 — 1CT 3 ) by the constraint from the exper- 
imental limit of the neutron EDM||. In our calculation, we always take \x to be real and thus 
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no CP-violating effects are induced by fi. The CP violating effects induced by Ad,u are greatly 
suppressed because they are proportional to the masses of the corresponding quarks which can 
be neglected compared with the squark mass parameters in £ sca i ar except that induced by A t 
which is associated with the top quark(see the form of squark mass matrices in (A. 8) and (A. 9) 
in Appendix A). Thus, arg(A t ) is the only new CP violating source in our calculation. 

After the interaction terms in the potential are diagonalized, the MSSM Lagrangian will be 
expressed by mass eigenstates instead of gauge eigenstates. The CP violating phases are then 
transferred to the gauge interaction vertices (see Appendix B for related Lagrangian pieces). This 
is reflected by mixing matrices in the interaction vertices. The mixing matrices Z + ,Z~,Zn which 
diagonalize charginos and neutralinos are real if [i is taken real. The mixing matrices Z t for the 
top squark is in general complex due to the complexity of A t . This implies that the CP violating 
effects comes from arg(A t ). The mixing matrices will be discussed in detail in the appendix A. 



3 Calculation 

We now discuss the CP violating effects in the process t — > bu T f, which corresponds to diagrams 
in Fig.l, within the MSSM framework. First denote the invariant mass of f and v T as y/q 2 , where 
q = p Vr + Pf . p Vr and p f are the four-momenta of v T and f. We calculated the CP asymmetry 
when q 2 > m^. This condition opens a new window so that x + an d X° cu t ma y gi ye an absorptive 
part to the amplitudes for loop diagrams as depicted in Fig. 2. 

Several points should be indicated at the moment. 

(1) . There should be a minus sign in front of the amplitudes for box diagrams relative to that for 
triangle diagrams. [] 

(2) . There are two kinds of box diagrams. One of them, Fig.2a and 2b, is fermion number 
nonconservative[[7| . 

(3) . We have assumed that the absorptive part of the amplitude for the loop diagrams are induced 
by X° X cut [ 12 1 . t or f can also be on shell and give new contribution to absorptive part under the 



condition m 4 ~ + m x o < m t or m x > rrif, respectively. We excluded these two cases for simplicity 
for the following reasons. Under our assumption about the stop mass, m 4 ~ + m x o < m t can be 
satisfied only in very narrow SUSY parameter region which is simply excluded in our calculation. 
m x > m f means the CP-odd effect will appear as \fcf > 160GeV (we take = lSOGeV^) which 
has too small branch ratio and can be ignored. 



1 This can be explained by that we can get the amplitude for Fig.2e from that for Fig. 2a by interchanging the 
field operators of \~ an( i ^ a t two ends of the f propagator and simultaneously changing f to W boson in the 
Wick contraction procedure. The interchange of two fermion operators gives the minus sign. An analogous case 
on QED can be found in EjJ. 
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Two quantities are defined to represent the CP-asymmetry effects, 

A t, e _ r-f 

A CP 



r + r 

JL t , e "g SP JLt, e a< i ~d^~ /yx 

L t = ( m x + m x o ) 2 , (8) 
L e = ( fOOGel/ ) 2 , (9) 

where - 1 and - 1 are differential widths of top quark and top anti-quark. A^p reflects the 
theoretical CP-odd effect appearing when the invariant mass yfq 1 is above the threshold. A e CP 
reflects the experimental CP-asymmetry effect when we measure the decay events with q 2 > L e , 
which is fixed to be (fOOGeV^) 2 in the work. 



We have only considered the contribution to the denominator in Eq. |7] from the tree-level 
diagrams which gives T = P. The nominator comes from the interference of the one-loop diagrams 
with the tree-level diagram, as 

Ar = r-f = A | M| 2 • phase space 

= ( |M| 2 - |M| 2 ) ■ phase space . (fO) 

The three body final state "phase space" is 

f f d 3 p b d 3 p f d 3 p Ur 



phase space 



2m t J (2n) 3 2E b (2n) 3 2E f (2tt) 3 2E Ut 

■(2-n)%p t -p b -p,-p VT ). (ff) 



The above expression multiplied by a S((p f +p Ul ) 2 — q 2 ) gives the phase space for fixed q 2 . M, M 
are the amplitudes for the process t — > bv T f and its CP conjugate process i — > bi7 T r, respectively. 

M can be expressed as 

M = aAt + 5^6*^, (12) 

i 

where the two terms come from tree-level and one-loop diagrams respectively, a, h % contain the 
CP-violating phases both from KM matrix element Vtb and from stop mixing matrix elements 
Zl 3 . A\ develops an absorptive part for q 2 beyond the threshold. M can be expressed as 

M = a*A 1 + Y,b i *Ai (13) 

i 

So, we have 

A|M| 2 = -4 J2 Im (a*b i )Im(AiA* 1 ). (14) 



ImA 2 is given by Cutkosky rule, 



]T VImAi = I/d* i(t- b X + X°) Mx + x° - fv T ) , (15) 

i 

where 

r <7 3 k H 3 V 

* = / wnE,, wm^ ' - Pb - h - k ' ] • (16) 

are the phase space of x°iX + as they are on shell, k, k' are the four-momenta of x° an d X + > 
respectively. After summing up all spins of external particles we get, e.g., the interference term 
of Fig. 2a and the tree-level graph of the form 

\ £ A|Af| 2 (a) = -2£/m(a*V) £ Im{A\A\) 
spin » spin 

-<K(g-*0 -raj ^ ^ (17) 



g 2 - m 2 w (2tt) 2 7 £ VVi 7 * 7 ^ P 2 

where we have introduced X*s to represent quantities like Im(a*b l ) in Eq. [14| arising from the 
SUSY couplings and the corresponding FjS to represent quantities like Im(A2A\) which are Lorentz 
invariant functions of the four-momenta of x° i b, f and u T . V\ and V2 are denominators of the 
two boson propagators in the loop graphs. Note that the SM phase from Vtb is cancelled in the 
interference of the tree diagram and the one-loop diagrams. 

We get the analytic expressions for |yj spin A|M| 2 by integrating the phase space of x° an d 
X + in q = system. The analytic results are given in appendix C. By expressing the formulae in 
Lorentz invariant form, we translate the formulae to the top quark rest system. In this system 
the final state three body phase space integration is implemented numerically. 

All the X l s for each graphs are proportional to £f (for detail expressions of X l s, see Appendix 
C), so, we can write 

Acp = & 1 • fcp, (18) 

where 

H = ImZt j *Z? = t±L Bin(20 t ) sin fa, (19) 

9t and fa are given in (A. 13) and (A. 9) of Appendix A. In terms of parameters in the top squark 
mass matrix we get that 

H = 1 ^ -7= 1 (20) 



From the expression of A in (A. 13) of Appendix A we can see that can be as large as | when 
the following conditions are satisfied at the same time, i.e., Lf = Rf, \i = and A t is purely 
imaginary. This is certainly difficult to reach. 
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4 Numerical results 



We now turn to our main numerical results. The calculation is based on the low energy MSSM 
scenario whose parameter freedom has been greatly reduced as described in section 2. Another 
simplifying assumption taken in our calculation is the universal relationship between the gaugino 
masses, i.e., m± = |m2tan 2 ^ 5 where 9w is the weak mixing angle flT3|] . We write the parameters 
m~ L , m~ R in the top squark mass matrix as 

m 2 iL = M 2 - cm 2 , m\ R = M 2 — 2cm 2 (21) 

where M is an arbitrary mass scale for scalar particles. 

Neglecting the masses of r and quarks except top quark we are left with ten SUSY parameters, 
i.e., [A, mi-, tan/5, c, M, \A t \, a.rg(A t ), mft ,mfR and mj. We take m^L = mfR — m u — 130GeV 
to which the results are insensitive, and always take \A t \ = M, c = 0.1 ~ 1. The other free 
SUSY parameters are restricted by the experimental limits on the masses of super particles |J 
and our assumption m t ~ > 140GeV. In particular, we take that, m x o , the mass of the lightest 
neutralino, is above 30GeV and m x + , the mass of the light chargino, is above 65GeV. Another 
limit for the parameter space is adopted for simplicity that we require m x o i + m x + l > lOOGeV^ 
and m x o i + m^ > m t . The SM parameters are taken as m t = 175GeV, \Vtb\ 2 — 1, a — 
m w = 80.33GeU , sin 6 W = 0.232 and m b = m T = 0. 

A consequence of the above scenario, especially that arg(/x)=0 and m b = is that Figs. 2c, 2d 
and 2f do not acquire any CP-violating phase. So only Figs. 2a, 2b and 2e are considered. It is 
found numerically that more than 90% of contributions to Acp come from the triangle diagram 
Fig.2e. Because of this, the results are not sensitive to the values of and mp. 

We have studied the CP asymmetry, A C p, as a function of SUSY parameters axg(A t ), /x, tan /3, 
m2 and 777^. In all the figures there are two curves for the same values of the fixed parameters 
of which the curve giving larger Acp represents A l C p while that giving a smaller value represents 

Si C p. 

In Fig. 3 we show the Acp as a function of arg(A t ). We can see that Acp is approximately a 
sine function of axg(A t ) as can be seen in Eq. [18] and Eq. 19. fcp defined in Eq. [L^ is plotted 
in Fig.4. We see that fcp is just like a parabola. As a result Aqp does not reach its maximum 
when A t is purely imaginary, rather, it is maximal at aig(A t ) w ±0.77r. fcp depends on axg(A t ) 
through the top squark mass (see Appendix A). 

In Figs. 5-7, we plotted A C p as a function of the Higgs mass parameter /i for tan/3 = 1.2, 5, 15 
respectively, for different values of 7772. The global feature of the three figures is that Aqp decrease 
dramatically as |/i| increase. Notice that, in the high tan /3 scenario Acp becomes quite insensitive 
to the sign of /1. Aqp is almost symmetric about \\x\ for tan/3 = 15. However, as tan/3 = 1.2 Acp 
is not small only for minus /x Another feature of the figures is that for fixed tan/3 and //, A e CP 
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decrease whereas A l CP increase as m 2 become larger. The reason is evidently that the threshold 
varies with m 2 . 

The dependence of A C p on tan/3 is plotted in Fig. 8 and Fig. 9, for several values of m 2 and for 
ji = —70GeV and ji = — 50GeV, respectively. An interesting feature of these two figures is that 
Aqp decreases as tan (3 increases whereas A e CP increases. This is because low tan f3 gives strong 
Yukawa coupling for top quark so that we get large A l CP . However, large tan (3 makes the threshold 
lower and thus elevates A e CP . So, A e CP prefers bigger tan (3 in contrast to the situation in the top 
quark two body decays. As tan (3 > 4, we can see from Fig. 8 that A l CP is almost insensitive to 
tan/5. 

The dependence of A C p on m 2 is plotted in Fig. 10, for several values of tan (3. We can see that 
the major part of the curve A e CP falls into the region between 0.1%-0.2%. For tan/3 = 2, A C p 
rises as m 2 increases whereas for large tan/5 A C p slightly drops as m 2 increases. 

Finally, we give the dependence of A CP on the top squark mass in Fig. 11 for /j, = — 50GeV. 
It is found that Acp depends essentially only on m^, not separately on M and c. A C p decreases 
with m ( -, just as what has been expected. When is around 140 GeV Aqp can reach up to 
0.5%. 

In summary, A e CP is around 0.1% level in the major part of the parameter space that we have 
discussed. In a very narrow region of the parameter space, A e CP can reach 0.5% level. 

5 Discussion 

We would like to point out two points in this section. 

(1) . The branch ratio for top quark decay drops rapidly when W boson is off shell. For 
yfq 1 > lOOGeV the branch ratio for bv T r final state is only about according to our calculation. 
However, as we have mentioned that the contributions to Aqp comes mainly from the triangle 
diagrams in Fig.2, Acp for different final states have the same sign and approximately the same 
size. So we can make a combining analysis for the data of the three body CP asymmetries for 
different three fermion final states (r, u T ), (//, v^), (e,v e ), and all allowed three quark decays . 
The total branch ratio for top quark three body decay for y/q* > lOOGeV can reach up to about 
Y^q. However, there are CP asymmetries of the order (^) 2 for three quark decays coming from 
gluon and gluino corrections, which may be of the same order as considered here. This needs 
to be studied further if one attempts to quantitatively compare the theory and experimental CP 
asymmetries for three quark decays. 

(2) . The total width of top and anti-top quark are equal T = f due to CPT theorem. The 
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following relation holds 

r(t - bff) - r(t - 6/7') = - ( r(t - b x \ + ) - f (t - frV) ) //; , (22) 

where //' represent all the W boson decay products mentioned above and the RHS of Eq. |22| 
denotes the CP asymmetry from the diagrams with //' as intermediate states. The relation can 
be easily seen from the Feynman diagrams representing these two processes in Fig.l, 2 and Fig. 12. 

The four interference terms contributing to the RHS of Eq. |22| between the two triangle dia- 
grams and the two tree-level diagrams in Fig. 12 correspond to the interference terms between the 
four box diagrams in Fig. 2 and the tree-level diagram in Fig.l. The other two interference terms 
of the fermion loop diagram with the two tree-level graphs in Fig. 12 correspond to those of the 
triangle diagrams in Fig. 2 with the tree-level diagrams in Fig.l. The relative minus sign on the 
RHS of Eq. ^ can be explained as following. We must add a minus sign to the amplitude for 
Fig.l2e due to the closed fermion loop. Such a minus sign is absent in the corresponding amplitude 
in Fig. 2 . On the contrary, there is a minus sign in front of the box diagram in Fig. 2. No such 
minus sign is present in Figl2.c, d. Thus, the CPT relation Eq. |22| manifests itself through the 
Cutkosky rule in Eq. [15]. 



6 Summary 

In this work we considered the CP asymmetries induced by MSSM new phase in the semi-leptonic 
three body decays of the top quark under such an assumption that the top squark is so heavy 
that no CP-odd effects are observable in top quark two body decays in one-loop level. In our 
calculation Acp can reach up to 0.55% in the most favorable case. Considering the small total 
branch ratio for yfq 1 > lOOGeV which is about ^5 f° r a ^ three body decays it is really hard to 
detect such small effects experimentally. 

However, several constraint in our calculations can be relaxed. For example, the gaugino mass 
parameters mi and mi can be complex and give new CP-violating sources if we do not make the 
universal assumption about the gauginos masses in GUT scale. Another possible CP source is 
from the Higgs mass parameter /1. According to recent studies on the neutron EDM, cancellation 
among different contributions can take place so that /x can have large imaginary part even when 
super particles are in O(100GeV) level [p^] . If this is the case, complex // can introduce additional 



CP asymmetries in top quark decays. And, if we relax the constraint on the to be slightly 
above lOOGeV, m% x and m x o cut can also give contribution to Aqp in three body decays. Another 
improvement may come from branch ratio enhancement. The direct experimental constraint on 
m x o + m x + is even below mw today 0. As yfq 1 decrease the branch ratio for top quark three 
body decay increase rapidly (e.g., the branch ratio for all three body decays is about ^ for 
x/q 1 > 90GeV). So, if m x o +m x + is not much heavier than m w , the A C p in the three body decay 
is hopefully detectable in the LHC which may be able to produce 10 7 tt pairs. 
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On the contrary, if m x + and m t ~ are both heavier than, e.g., 140GeV, all the windows for 
CP asymmetries induced by super particles will be shut up. CP-odd effects in top quark decays 
induced by MSSM particles can only exist beyond one-loop order. This will be beyond the 
experiment ability in the near future. 



Appendix A 



There are different conventions for super particles mass matrices adopted in literature so that 
it is very easy to make sign errors. To avoid sign errors we rederived the MSSM Lagrangian. 
For most part we adopted the same conventions described by J. Rosiek||. The original MSSM 
Lagrangian is given by Eq. [I] — Eq. || of which all fields are gauge eigenstates. We should point 
out that the eu = —1 convention will give an extra minus sign to parameter fi compared to 
those adopted the convention ei2 = 1 (if universal relation for gauginos is taken only relative sign 
between the gauginos and \i is significant). To get the physical spectrum of particles one should 
carry out the standard procedure of gauge symmetry breaking. After SSB super-particles will 
mix and form different mass eigenstates. 

The charged Higgsinos and charged winos mix and give two mass eigenstates named charginos. 
The mass matrix of charginos is 



m 2 \/2Mw sin f3 

V2M W cos p fi 



The mixing matrices satisfy 



and is defined by 



[Z-) T M X Z + = diag (m Xl , m 




¥2 



(A.l) 

(A.2) 
(A3) 



-i A" 

ip 1 



H 2 
2 




(A.4) 



In the above equations A ± = ^{^w-fi^w) where A^ are the first and second components of wino 
as in Eq. [5[ is the first (or up) Fermion component of the second Higgs super field doublet. 
The fields on the left hand side of the above equations are gauge eigenstates and the fields on the 
right hand side are mass eigenstates. The four-component Dirac spinor charginos are defined by 



<Pi 



The mass term which will appear in the final form of Lagrangian is — m Xi XiXi 
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The third component of wino, photino and neutral Higgsinos combine to give four Majarana 
neutralinos. The mass matrix for neutralinos is 



M x o = 



mi 

1712 

—Mz cos (3 sin 9w Mz cos (3 cos 9w 

M z sin P sin 9 W —M z sva.f3 cos 9 w 



which is diagonalized by 



— M z cos (3 sin 9w 
Mz cos (3 cos 9w 


-n 



M z sin (3 sin 9w 
—Mz sin (3 cos 9w 
-fi 




(A.5) 



Zn is defined by 



ZjfM x oZ N = diag (m x o, 







(A.6) 



1>l 



V ^ 2 



Hi 

2 

#2 



= z 



(<p\\ 

v y?4 y 



(A.7) 



where A B is photino. All fields on the left hand side of the above equation are gauge eigenstate 
and those on the right hand side are mass eigenstates. The four-component Dirac spinor form 



for neutralino is Xi — 



The mass term in the final form of Lagrangian of neutralino is 



Ignoring generation mixing, the mass eigenstates of squarks are obtained by mixing the left- 
handed and right-handed eigenstates of squarks. Its mass matrix is 



M) 



Rf 

c f 



J S 



rn 



m 



Cf 
C) R f 

2 
/ 
2 



(A.8) 



cos 2(3(T 3f - Q f sin 2 9 W )M 2 + mj L 
f -r cos 2(3Q f sin 2 9y/M\ + m 2 ^ 
-m f (r f fx + A)) = \C s \j*s 



(A.9) 



where T 3 j is | for up squark and — | for down squark. Qf is the charge of the sparticle and rf 
is cot (3 for up squark and tan/3 for down squark. m 2 ? , m| , Af are the corresponding diagonal 

JL JR 

elements of wig, rrifj^mjj) and Ay (Ad) in £ so f t , respectively. The mixing matrix satisfies 



ZfMjZf = diag ( 
They are defined, e.g., for up squark, by 



i i 

171: ,171; 



u 



R 




(A.10) 



(A.ll) 
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where Q U L is the up component of left-handed up squark doublet, U R is the right-handed up squark. 
Ui t 2 are the two mass eigenstates of up squark. The conjugate of U R in Eq. |A.11| comes from that 
we adopt the charge conjugate of left hand Fermion to represent its right hand component in the 
original MSSM Lagrangian. The final form of squark mass term in Lagrangian is —mfj.Ui Ui. In 
particular, the mixing matrix for top squark is given by 



cos 9 t e i(t>t/2 - ■ 



sin 6 t e~^ t/2 cos t e -i **/ 2 



(A.12) 



where 4> t is defined in (A. 9) and 



tan 9 t 
A 



2\Ct 



Li- 
ft 



R, 



Rtf + MCt 



(A.13) 



The formulae for down squark are similar to up squark. Note the definition of mixing matrix for 
down squark given by || is the complex conjugate of the mixing matrix given here. 



Appendix B 



In this appendix we list the relevant pieces of the SUSY Lagrangian in terms of the mass 
eigenstates ||. 





= gt*x^P L + B^P R ]t + H.C. 


(B.l 


£-bt x 


= gt i b[C^P L + D^P R ]V t lx- + H.C. 


(B.2 


£bb X ° 


= gb*x° j [E ij P L + F IJ P R }b + H.C. 


(B.3 


^tb X 


= gb*Xj [G ij P L + H*P R )V£t + H.C. 


(B.4 




= gf*x° 3 [ MlJp L + N l3 P R }r + H.C. 


(B.5 


r - 

t~>VT X 


= gU l it*x 3 P L v T + H.C. 


(B.6 


r - 

'-'TVX 


= -gv* T x~2[z+p L + rz^*p R ) + H.C. 


(B.7 


£vux° 


= gW i i>W i P L v T + H.C. 


(B.8 


C xx o w 


= 9X^[O ij PL + V ij P R ]x W+ + H.C. 


(B.9 
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where 



Z li *az% sm6 + Z 2 J cos 9) - u l Zf*Z% 



\/2 cos C 

2y / 2sin6> y2i* ylj* , ,t yli* yij* 
3 cos 9 Z < Z Af t Z W 

-d h ZfZ^ 
-ZfZt* +u t ZfZ 23 

-^- Q Zl%Z^ sin - Z# cos 0) + d% 2i i$ 



Vising y2i yli* , ifc yli ySj* 
Q„„„a ^/V "T U Zj h Zj n 




d h ZfZ 2d 



(B.10) 



-jJ—^iZ^ sm8 + Z 2 J cos6) + FZ?Z 3 J 

-Vising y2i yjj* _|_ \ T 2j U Z^* 



1 y 4 iy- 

Z2i* y 1 1 y3i* y— 
at + 7f ^2 

it*, c? b and Z r are the Yukawa coupling parameters for top, bottom and r respectively. They are 
given by 

«' = -r"^ > 

V2mvi/ sin p 

d * = ZOO* 

\/2mw cos (3 

r = ~ mr - ■ (B.n) 

\/ lm w cos p 



The vertex in Eq. [B.9| can be expressed as 

£>xx°w = -9 X y[O lJ *P R + V^Pl^W- + H.C. (B.12) 
which can be used to read off the Feynman rule directly for the Fig.2e. 

Appendix C 

In this section we give the formulae for the quantity ^yjA|M| 2 in Eq. [TJ] for Fig. 2c and 
2e in detail. In the following formulae the X i s are not separated out explicitly as in Eq. |D] for 
convenience. All the formulae and variables are given in the q = system. The imaginary part of 
the amplitude for Fig2.c is 
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u{p u )T 9ux o{k' + m x o)T Ux ou{p t )u{p b )T ibx {i m x )T DTX v{p T ) 
l(p t - k) 2 - mf] [{p u - kf -ml] 

where k = (E, k) is the four-momenta of x° and the Ts are interaction vertex factors. The quantity 
|E spm A|M| 2 forFig.2cis 

1^ ,,,,,,, , o 6 \k\ { 1 , m r i . r*. . 2p6 cos a 

-£A|M| 2 (c) = y ' ' — (2r + L^|fclcosa) 1 



2^ 1 lw q 2 - mly 47TyV I A' a ^' ' y A'S'a 

a(S + LS'/4) + T + iL£?' \k\ cosa, 1 + a 

log- 



A' a 2 to 1 - a 

1 / O Z 1 1 \\ -, 1-6 



p. 1 , l-a&cosa + v^l ,„ „ x 

+ 7i)^77 lo g; ; 777 ( C - 2 ) 



A'B' a a 2 ' \[K I- abcosa- 



where 



o v2 I ; 1 2 1 -m2 

P = -2X a Pv ■ PT\k\ \Pt\ 

1 o — * 

Q = (A' a m x om t -2A' a (p b -p J/ + A;-5))p I ,-p T |A;||p t | 

+4* a V • Pr^l^l IpII - 2tX t £ - b't 

£ = 2A' a 1 m x om t p fc • • p r + {2X 2 (p b ■ p u + k ■ q) - X^m x om t )p u ■ p T EE t 
+Xa m x om tPv -Prk-q- 2X 2 p„ ■ p T (EE t ) 2 

+2S'E' Ur E - le' Vt e 2 + b's + LB' /A 

-X 2 (Pb ■ VvVt -Pr-Pb- PtPu -Pr+Pf PuPb ' Pr) 

K = a 2 + b 2 -a 2 'b 2 sin 2 a- 2ab cos a (C.3) 

In above formulae, E' Vr is the energy of p v in the q = system, a is the angle between the 
three-momenta of top quark and u T in the q = system, q 2 is the invariant mass square of the 
final state lepton pair. Other quantities are defined as: 



A 2,2 2 

A = m t + m x o — rrif , 

B = m x o — to~, 

A' = A — 2EE t , 

B' = B-2EE' Ut , 



a 



b = 



A' ' 
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L = X aPfPb, 

S = 2 X ! m x om tPb ■ q + Xjpb ■ PvPt ■ q + 

*a(Pb -Pr-Pf Pv)EE t + X 2 p t ■ p v k ■ q, 
T = X 2 (p b -p T -p t -p„)\k\\p t \, 

S ' = 2 X l m x om t(Pb -Pr-Pb- Pv) + X%p b ■ p v p t ■ q + 

X a(Pb -Pr+Pf Pv)EE t - X 2 p t ■ p u k ■ q, 
T' = X 2 (p b ■ Pr + Pt ■ Pv)\k\\Pt\ (C.4) 
The non-Lorentz invariant four- vector component are expressed as 



E t 

\ft\ 
E 



m 2 + q 2 

i i 
m t — q 

m 2 + q 2 
2m t 



\k\ = JE 2 -m 



0, 



m\ + q 2 - 4p t ■ p u 

cos a = k t, (^-S) 

mf — q 2 

Xl' 2 are given by 

Xl = -ImB^D^Z+W 1 * , X 2 = -ImA ll D il Z^W l * (C.6) 

The elements of the mixing matrices in the above equation are taken according to that only m x o 
and m + are considered in our calculations. 

The corresponding formulae for Fig2.e are 

- 4 - d 3 k 



ImA ( e ) = — ^ / ^5((q - k) 2 - m 2 ) 

u(p b )r fbx (U' + m x )-t tx Y xx0w V itx0 u{p t )u{p v )^ lx P L v{p T ) 



(q 2 - m^lipt - k) 2 - m~] 
The quantity \ £ spi „ A|M| 2 for fig.2e is 



(C.7) 



2^ 1 1 ' (q 2 - m 2 w ) 2 2V2n 

U 2( y_^ ) + (Z _r + 4 )log i±^\ (c.8, 
Lai a a a z 1 — a) 
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where 

X = 2Xl[Ef Vr \P t \ cos a(p b -p T -p t -p„)-\P t \ 2 p v -p T + 

~E' v \os 2 apfp b ]\k\ 2 

Y = VI- S\P t \\k\+TE' UT cos a 
Z = £ + SEE t + T'EE' Ut + H 
S = (2X^p b ■ p u - X^m x om t )p u ■ p T 

T = 2X 2 m x m t p b ■ p v - X*m x om t p b ■ p t - 2Xj?q ■ p t p b ■ p u 
T' = 2X 2 m x m t p b ■ p u - X^m x om t (p b ■ p u - p b ■ p T ) - 

2Xlp b ■ p v {p t -Pu-Pt-Pr) 
n = 2X^ {-p b ■ p T EEl^{\P t \+E t cos a) + Pt -p u EE' T \k\(E t cos a -\P t \) 

+(2EE t - k ■ q)\P t \\k\Pu ■ p T ~ k ■ qE' Vr \k\ cosa Q ■ p v } 
£ = 2X e 3 {E 2 E t El T {p b -p T +pt-Pu) -EE' T k-qp t ■ p u - 

(EE t - k ■ q)EE tPu ■ p T - E' V \E 2 - ^\k\ 2 )p t ■ p b } 

H = X^m x om t p u • p T k ■ q + 2(X^m x om x p t ■ p T + X^m x om t p v ■ p T )Pb • Pv 

-X^m x o (p b ■ p v p t -Pr-Pb- PtPv -Pr+Pf PuPb ■ Pr) (C.9) 

The parameters X % e s are given by 



Xl = ImA il D il (-V 11 *), X 2 = ImB ll D ll {-V 11 * 
Xl = ImA* 1 D il {-O n ), Xt = ImB il D il \-0 



(CIO) 



From the above expressions and the expressions for A, D in Eq. p.!0| we can see that the Xs are 
proportional to Q = Im{Z\ l * Z 2% ) as pointed out in section 3. 
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Figure Captions 



FIG. 1 The tree-level Feynman diagram for the process t — > bu T f. 

FIG. 2 The SUSY induced CP-violating one-loop diagrams for the process t — > hv T f. 

FIG. 3 The CP asymmetry Aqp is plotted as a function of axg(A t ) for tan (3 = 1.2, 777 2 = 150GeV, 
H = -AOGeV, M = 200GeV, c=0.2. When m- tl = 150GeV, A CP reaches its maximum. 

FIG. 4 The quantity fcp defined in Eq. [18] is plotted as function of axg(A t ). All the parameters 
are the same as that of Fig. 3. 

FIG. 5 The CP asymmetry Acp is plotted as a function of SUSY parameter ft, for several values 
of m 2 , for tan (3 = 1.2, M = WOGeV, c=0.15, axg(A t ) = 0.5tt. 

FIG. 6 The CP asymmetry Acp is plotted as a function of SUSY parameter ft, for several values 
of 777.2, tan/? = 5. All the other parameters are the same as that of Fig. 3. 

FIG. 7 The CP asymmetry Acp is plotted as a function of SUSY parameter /j, for several values 
of 777.2, tan/3 = 15. All the other parameters are the same as that of Fig.3. 

FIG. 8 The CP asymmetry Acp is plotted as a function of tan/3, for several values of 777-2, 
H = -70GeV, M = WOGeV, c=0.15, arg(A t ) = 0.5tt. 

FIG. 9 The CP asymmetry Acp is plotted as a function of tan/?, for several values of 7772, 
fi = -50GeV, M = 160GeU, c=0.15, arg(A t ) = 0.5tt. 

FIG. 10 The CP asymmetry Acp is plotted as a function of 777 2 for M = 160GeV, c=0.15, 
arg(^4t) = 0.57T, fi = — 50GeV for tan (3 = 2 and ft = —60Gev for tan f3 = 5, 10. 

FIG. 11 The CP asymmetry Aqp is plotted as a function of the mass of the light top squark, for 
different values of 7772, tan/3 = 2.5, fi = — 50GeV, a.ig(A t ) = 0.57T. 

FIG. 12 The tree-level and one-loop Feynman diagrams for the process t — > b\°X + ■ 
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